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Abstract 

Using the Landauer formula and a random matrix model, we inves- 
tigate the autocorrelation function of the conductance versus magnetic 
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field strength for ballistic electron transport through few-channel mi- 
crostructures with the shape of a classically chaotic billiard coupled 
to ideal leads. This function depends on the total number M of chan- 
nels and the parameter t which measures the difference in magnetic 
field strengths. Using the supersymmetry technique, we calculate for 
any value of M the leading terms of the asymptotic expansion for 
small t. We pay particular attention to the evaluation of the bound- 
ary terms. For small values of M, we supplement this analytical study 
by a numerical simulation. We compare our results with the squared 
Lorentzian suggested by semiclassical theory and valid for large M. 
For small M, we present evidence for non-analytic behavior of the 
autocorrelation function at t = 0. 

PACS numbers: 72.20.My,05.45.+b,72.20.Dp 
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Recent advances in semiconductor technology have made it possible to 
study experimentally the transport of ballistic electrons through microstruc- 
tures [I], with the shape of a classically chaotic billiard. Two aspects of the 
dependence of the conductance G(B) = (e 2 /h)g(B) on an external magnetic 
field B have received |3], |J] particular attention: The suppression of weak 
localization at small values of B, and the form of the conductance autocor- 
relation function C(AB) = 5g(B)5g(B + AB). Here, 5g(B) = g(B) - ~g(B) 
is the difference between the dimensionless conductance g(B) and its mean 
value. The bars indicate an average which experimentally is taken over the 
Fermi energy or the applied gate voltage. 

In a previous paper |J , we have calculated the generic form of the weak 
localization peak. In the present paper, we address the generic form of the 
conductance autocorrelation function C. For ideal coupling between leads 
and microstructure, this function depends on the number M of channels in 
the leads, and on t, a measure of AB to be defined below. 

Several authors have calculated the autocorrelation function C in an ap- 
proximate way. In the semiclassical approximation ||, C was found to have 
the form of a squared Lorentzian. The semiclassical limit is expected to ap- 
ply only when M 3> 1. Moreover, the non-diagonal terms in the Gutzwiller 
trace formula are neglected. Attempts to determine C from a random matrix 
model for the scattering matrix with the help of a Brownian motion model 
[0, H have not been successful. Use of a random matrix model and the su- 
persymmetry technique || [10[ has confirmed the squared Lorentzian shape 
of C. However, these were also asymptotic calculations valid in the limit 
M»l and thus kin to the semiclassical approximation. 

In the actual experiments || ||, the number of channels per lead was quite 
small. Thus, it is not clear whether the theoretical results just mentioned ap- 
ply to this situation. Moreover, a theoretical study of the full autocorrelation 
function is of considerable interest in its own right. Indeed, C is a four-point 
function, and to the best of our knowledge it has not been possible so far 
to calculate this kind of function without approximation analytically, using 
a random matrix model as starting point. The work reported in this paper 
does not solve this problem fully either. However, we succeed in calculating 
analytically the leading terms in an asymptotic expansion of C for small t. 
We believe that in so doing, we display interesting and novel aspects of the 
supersymmetry technique, especially in relation to the boundary terms. In 
addition, we are led to question whether C is analytic in t near t = 0. 
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Starting point is the random matrix model for the Hamiltonians of the 
chaotic microstructure for two values of the external field. Our approach 
is very close to that taken in Refs. R 110], and also to our work on weak 



localization ||, except that we now have to deal with a four-point function. 
We use the Landauer formula and Efetov's supersymmetry technique |T]J in 
the form of Ref. [0. 

The calculation of the leading terms of the asymptotic expansion of the 
autocorrelation function for small t involves a generalization of the usual po- 
lar coordinates used for two-point functions. The transformation to these 
variables necessarily introduces the boundary or Efetov-Wegner terms into 
the formalism. We devote particular attention to these terms. To the best of 
our knowledge, integrals in polar coordinates over graded vectors have first 



been treated by Parisi and Sourlas [13], and over graded matrices by Efetov 



|TT| ] and Wegner, see Ref. [F2 . Our central result given in Eq. ( [71] ) below 



has first been derived in a different way by Efetov [ITJ and by Zirnbauer [[TJ 
The derivation which we employ makes use of the method of boundary func- 
tions suggested by Berezin ||15||. For a more general discussion of coordinate 



transformations in graded integrals we refer to the papers by Rothstein [16 



and by Zirnbauer and Haldane [T7\. 



As a complement to the analytical work which yields only the leading 
terms of the asymptotic expansion, we simulate the correlation function nu- 
merically. We focus attention on small values of the channel number M. 

In Section [l], we formulate our approach. The integral theorem used in 
calculating the asymptotic expansion, is derived in Section |2|. The analytical 
calculation is discussed in Section ||, our numerical results are presented 
in Section |j. Section |^ contains the conclusions. Additional mathematical 
details are given in two Appendices. 



1 Formulation of the Problem 

Much of the development in this section is similar to that of ref. 0] where 
further details may be found. We try to keep this presentation both brief 
and self-contained. 
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1.1 Physical Assumptions 

We thus consider ballistic electron transport through a two-dimensional mi- 
crostructure with the shape of a classically chaotic billiard in the presence of 
an external magnetic field B. We wish to determine the generic features of 
the autocorrelation function 

6g(BW)6g(BW) = g(BW)g(BW) -gJWtfgjBWj. (1) 

Here, g(B) is the dimensionless conductance, and the bar denotes an average. 
Experimentally, this average is typically taken over the applied gate voltage 
[[|, while theoretically we take it over a suitably chosen random matrix model. 
We assume that the Hamiltonians of the closed billiard exposed to the 
two fields B^ r \ r = 1, 2, are described by two ensembles of Hermitean random 
matrices HffJ of dimension iV which have the form 

H$ = H$-^H<ZP, H$ = H$ + ] J±HW. (2) 

The matrices Hj£p with C = I, II belong to two independent (uncorrelated) 
Gaussian unitary ensembles (GUE's) with first and second moments given 
by ^ 

Hf2 = 0, Hf£ hJ^J, = —SfMs'Sii'w (3) 

The parameter A fixes the mean level spacing d of each ensemble and must 
be adjusted to the local mean level spacing at the Fermi energy Ep of the 
billiard. In the limit N —>■ oo always considered in this paper, the average 
level density d~ l takes the shape of Wigner's semicircle, and the Fermi energy 
Ep = is chosen to lie at its center where d = ir\/N. As shown in ref. [0], 
the parameter t is related to the area A of the billiard by \/i = k\B^ — 
B^\A/(2<p Q ) where k is a numerical constant of order unity and where 0o = 
hc/e is the elementary flux quantum. 

For this model to be valid, the fields B^) must obey the following two 
conditions, (i) B^ must be strong enough so that for both r = 1,2, the 
crossover transition from orthogonal to unitary symmetry is complete. In 
practice this requires fields stronger than a few millitesla, cf. Ref. ||. (ii) 
The fields must be weak enough for a GUE description to be valid. In other 
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words, the classical cyclotron radius must be large compared to the diameter 
of the billiard. 

The billiard is attached to two quasi one-dimensional leads. In each lead, 
there are M/2 transverse modes which define the incoming and outgoing 
channels. These channels are labelled a = 1, . . . , M/2 in lead one and b = 
M/2 + 1, . . . , M in lead two. With \fi > denoting a complete orthonormal set 
of basis states in the billiard, and c the channels in either lead, we denote the 
coupling matrix elements between billiard and leads by W CfM (E). We assume 
that these matrix elements obey time-reversal symmetry, W Cfl (E) = W*^(E) 
so that the effect of the magnetic field is confined entirely to the interior 
of the billiard. The characteristic dependence on energy E of the elements 
W ctl {E) is typically very slow on the scale of the mean level spacing d, and 
is therefore neglected. 

We write g(B^) = gW and use for g^'\ r — 1,2, the Landauer formula 
for zero temperature, 



M/2 M 





a=l b=M/2+l 



The elements of the scattering matrix S t 



are given by 



Sl r J = 5 cd -2 i n^W c ,[D^X 1 u W du 



(5) 




(6) 



c 



We must put E = E F and recall that E F = 0. For a ^ b, we have 





where the trace runs over the level index /i, and where 




(8) 



For later use, we also define 



(9) 



c c 
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Because of the unitary symmetry of our matrix ensemble, the coupling 
matrix elements W CfJ- appear in all ensemble averages only in the form of 
the unitary invariants X cd = 7rX)„ W^W^. As in ref. |5[], we assume that 
X cd = X c 5 c d is diagonal in the channel indices. The parameters X c appear 
in the final expressions only via the "transmission coefficients" or "sticking 
probabilities" T c given by 

T C = 4AX C (A + X C )- 2 . (10) 

With < T c < 1, the transmission coefficients measure the strength of the 
coupling between billiard and leads. Motivated by the absence of barriers 
between billiard and leads, we put T c = 1 for all c (maximal coupling). Our 
model then depends only on M and on — B^\, i.e., on two parameters 
which are fixed by experiment. All dependence on Ep and on A disappears 
(the parameter A appears only in the transmission coefficients T c ). 

The ensemble average of g^ is obviously the same as for the case of a 
single GUE, = M/2 for r = 1,2. The calculation of the autocorrelation 
function, Eq. ([!]), therefore reduces to that of the ensemble average g^g( 2 \ 
To find g(^g( 2 \ we use both analytical and numerical methods. Our an- 
alytical study of the autocorrelation function presented in this and in the 
following three sections is based on Efetov's supersymmetry method WT 



in 



the version of ref. [12|. We refer to the latter paper as VWZ. The numerical 



calculation of the autocorrelation function is discussed in Section f|. 



1.2 Supersymmetry Approach 

We deal with a four-point function. This necessitates some changes in the 
formulation of VWZ. 

The generating function Z(J) is given by the graded Gaussian integral 

Z{J) = f X>[$] exp ( i £ < $ Qrp , [(D + J)$] Qrp > J (11) 

\ arp ) 

where we use the notation 

<F,G>=J2F(v)G(ri. (12) 
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The 8N dimensional graded vector $ has components $Q, rp (/j). The su- 
persymmetry index a distinguishes ordinary complex (commuting) integra- 
tion variables (a = 0) and Grassmannn (anticommuting) variables (a = 1). 
Later, the two values and 1 of the supersymmetry index a will also be 



denoted by b and /, respectively. As in subsection [TTT], the index r = 1,2 
refers to the two fields. The index p — 1,2 refers to the retarded and the 
advanced propagator as usual. Finally, /x = 1, . . . , N. The adjoint graded 
vector is defined by = &([i)s, with 

s , , , - A ,A ,A , (11) 

•->arp,a'r'p' \ J u aa' u rr' u pp' ■ \ ± ' J J 

The 8N dimensional graded matrix D is the graded inverse propagator, 

D = El 8 xl N -H + iLxQ, (14) 

with 

H = 1 8 X ff« - T 3 x . (15) 

The 8N dimensional graded matrices L and r 3 are given by 

Larp,a'r'p' ( ) ^aa'^rr'^pp' i [^~3\arp,a'r'p' ( ) ^ ^aa'^rr'^pp' • ^) 

The 8iV dimensional graded matrix J = J(e), a function of M 2 ordinary 
variables e^, has the form 

J(e) =J2l {rp) x J (rp) (e) , (17) 

where 



Ia'r^p',a"r"p" ~ ( — ) < - 1+Q ^ 5 rr /<5 p p/<5 a / a //<5 r / r »5py/ , (18) 

4v'(e)=^EWVWV47 } - (19) 



a'r'p ' ,a"r"p" 

rM/.-i ^^ir ,ir. ..,- ,r '" 

aft 

The integration measure is = flarp^ d$ Q rp(Ai)d$* rp (/i). Performing the 
Gaussian integration, we obtain 



Z(J) =Detg-\D + J) (20) 



where the symbol Detg stands for the graded determinant. Eqs. ([14]) to (|19|) 
show that the matrix D + J consists of four diagonal blocks corresponding to 
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the two different magnetic fields and to retarded and advanced propagators, 
respectively. The function Z(J) factorizes accordingly. Differentiating the 
resulting expression with respect to the source parameters e and averaging 
over the ensemble yields 

d 4 



\Sab Pl^o'l/I 2 - (11) (12) „ (21) „ (22)^(^)l e=0 • { " 2[) 



For later developments, we note that Z(J) possesses an important symme- 
try. We consider transformations T : — > T$(/x) which leave the bilinear 
form invariant. Such transformations satisfy the condition 

T- 1 = sTh . (22) 

The set of 8 dimensional matrices T forms a graded unitary group G = 
U(2,2/4) with compact and non-compact bosonic subgroups. For t = Q = 
J = 0, the generating function Z(J) is invariant under this group. 

1.3 Saddle— Point Approximation 

We average Z(J) over the ensemble, apply the Hubbard-Stratonovich trans- 
formation as usual to reduce the number of independent integration variables, 
and use the limit iV — > oo to define a saddle-point approximation. 
In Eq. ([TI]), we have to average the term 

exp I -i < ®arp, (H$) arp > I . (23) 

\ arp / 

Using Eq. (|2|), we find that the average is the product of two terms, 



ex P (-* E < $ <*n» (!§ x HW^) arp = exp (-^ tr gS 2 ) , ( 24 ) 



2 



\2l 



exp < ®<*r P , (r 3 x ^_ #(")$) arp > j = exp ^-__trg[(Sr 3 ) 

(25) 

where 

Sarp,a'r'p' *^ ^orp; ^ct'r'p' ^* ^ ] ^otp(aO ^a'r'p' (aO ■ (26) 
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Collecting results, we have with E = Ep = 0, 

Z(jj = J exp (iC(S) +iJ2< ®*r P , [{iL xn+ J)$} arp > j , (27) 



arp 



where 

iC(5) = -^(trg^ 2 + ^trg[(Sr 3 ) 2 ]) . (28) 

The following steps are quite standard and are only sketched here. After the 
Hubbard-Stratonovich transformation and integration over the variables $, 
Z(J) is written as an integral over the eight-dimensional graded matrices a, 

W)=J V[<t\ exp (-y (trg a 2 + l t rg[(ar 3 ) 2 ])) x 

exp (-Trg In [-X(a + ^r 3 ar 3 ) x 1 N + iL x Q + J]) . (29) 

The symbol trg (Trg) stands for the graded trace over matrices of dimension 
8 (SN), respectively. The saddle-point condition yields for the saddle-point 
manifold a G the result 



-iQ, Q = T~ l LT . (30) 



The matrices T obey Eq. (p2|) and belong to the coset space G/P defined 
by G/P = U(2,2/4)/U(2/2) x U(2/2), where P = U(2/2) x U(2/2) is the 
subgroup of matrices in G which commute with L. After integration over 
the massive modes and taking iV — > oo, Z(J) takes the form 

Z[T) = J V[Q) exp (iC eS (Q) + tC SO mce(Q, J)) , (31) 



where 

C eS(Q) = \ tr S t(^ r 3) 2 ] - Trg In [i\(Q + ^t 3 Qt 3 ) xl N + tLxQ] (32) 



and 

^source (Q, J) 



-Trg In 



Is x 1 A ( i\(Q + -^t 3 Qt 3 ) xl N + iLxQ\ J 



.(33) 
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Summing over the level index /i, introducing the transmission coeffients T c 
of Eq. ( |10|) and recalling that T c = 1 for all c, we obtain 

WW ] = J V[Q) detg(l + QL)- M R(Q) exp (-± trg[(Qr 3 ) 2 ]) , (34) 

where the source term R(Q) has the form 

R(Q) = 1m 2 tTg[G^G^]tig[G^G^} 
+^M 4 trg^^G^ltrg^ 21 ^ 22 )] 

+^M 3 trglG^G^G^G^ +G {11 ^G^G^G^] , (35) 



with 



G (rp) = G (rp)(Q) = + g^)-i/M . (36) 



The autocorrelation function is a function of two parameters, the channel 
number M and the parameter t for the magnetic field, 



C{t, M) = SgWSg® = gWgW - gW g <A , (37) 



with g^g( 2 > given by the graded integral (34). The graded matrices Q can 
be parametrized in terms of 32 real variables, half of them commuting, the 
others, anticommuting. In calculations of the average two-point function, 
one typically deals with a total of 16 integration variables. Except for this 
increase in the number of variables and for the form of the source terms 
(which are, of course, specific to our problem), the form of our result in 
Eqs. (j3~7D and (0) is quite standard. In spite of this similarity, the increase 
in the number of variables renders a full analytical evaluation of the graded 
integral (0) very difficult. As pointed out in the Introduction, our analytical 
work is restricted to the region of small t. We evaluate the leading terms 
of the asymptotic expansion of C (t, M) in powers of t. Progress in this 
calculation depends crucially on the proper choice of the 32 variables used in 
parametrizing the matrices Q. 
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2 Parametrization of the Manifold 



The integral ([34]) extends over the manifold of Q = T _1 LT, where T stands 
for the graded matrices of the coset space U(2, 2/4)/U(2/2) x U(2/2). Expo- 
nentiating the coset generators shows that the matrices T can be parametrized 
by two four- dimensional graded matrices £12 and £21, 



£21 VI + ^21^12 



(38) 



related to each other by 



£21 = kt{ 2 , (39) 

with k = diag(l, 1, — 1, — 1). (In presenting the matrices explicitly, we as- 
sume that the row labels of matrices of dimensions 8, 4, and 2 are (p,a,r), 
(a,r), and by r, respectively, and that the indices follow in lexicographical 
order. The matrices of dimensions 8 and 4 are presented in block form.) 
The 32 matrix elements of £12 and £21 represent the Cartesian coordinates 
of T. The range of the commuting coordinates is defined by the condition 
that the ordinary part of the matrix 1 + £12^21 be positive definite (non- 
compact parametrization in Boson-Boson blocks, compact parametrization 
in Fermion-Fermion blocks). Following the arguments of section 5.5 of VWZ, 
we find that the integration measure T>[Q] is identical with the invariant in- 
tegration measure for integration over the coset space. In Cartesian coordi- 
nates, this measure takes the form 

d//(T(W 21 )) = II d[t 12 \ar, a 'r'd[t 12 }* ar>alr , , (40) 
ar,a'r' 

where for ordinary complex z = x + iy, our dzdz* = 2idxdy. We write 
d/i(T) for d/i(T(ti2, £21))- Despite the simplicity of the measure, the Cartesian 
coordinates are not well suited for evaluating our integrals [0. Therefore, 
we follow Efetov and change to polar coordinates. 



2.1 Polar Coordinates 

As in VWZ, we express the matrices £12^21 and £21^12 in terms of their common 
eigenvalue matrix e and of the eigenvector matrices u\ and 112, 

£l2^21 = UiEU^ 1 , t 2 \t\2 = U26U2 1 . (41) 
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The unitarity relation ( |39| ) requires that both diagonalizing matrices u\ and 
u 2 belong to graded unitary groups of type U(2/2), 

u^ 1 = u\, «2 1 = ku\k . (42) 

To specify the matrices e, U\ and u 2 unambiguously, we order the eigenvalues 
e Q i and e a2 by the magnitudes of their ordinary parts, ord|ebi| < ordje^l, 
ord|e/i| < ord|e/2|- Moreover, we assume that the diagonalizing matrices U\, 
u 2 are elements of the cosets U(2/2)/U(l) x U(l) x U(l) x U(l), given by 

Ml = UixUxpUij, u 2 = u 2x u w u 2l . (43) 

Here 
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(44) 

are exponentials of coset generators specified by 8 commuting parameters 
£ pQ , £p Q , and by 16 anticommuting parameters j3 pq , j3 pq , 7 pr and 7* r , with 
a — b,f, and p, q, r = 1,2. We parametrize the matrices u px by their off- 
diagonal matrix elements x pa = e* arg ^ pa sin \£ pa \. We note that the matrices 
u{ x t\ 2 u 2 and u 2 t 2 \U\ commute with the eigenvalue matrix e, denote these 
two diagonal matrices by A and A, respectively, and write the matrices t\ 2 
and £ 2 i as products, 

ti2 = Ui\u 2 1 , t 2 \ = m 2 A% 1 . (45) 
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According to Eqs. ( p9[ ) and (f4"T|), the matrices A and A are related to each 
other and to the eigenvalue matrix e by 

A = kX\ e = AA. (46) 

Substituting Eq.§|) in Eq.fl38|) yields 

T = UAU-\ U = U x UbU~, 




with 

U w = diag(u lw ,u 2w ), w = x,(3,j, A = diag(A w , X b2 , X/i, A/ 2 ) . (48) 

These equations define the desired polar parametrization of T. Out of 32 
polar coordinates, 24 reside in the matrix U (parameters x pa , f3 pq , 7 pr , with 
a = b, f and p,q,r = 1, 2, and their complex conjugates), 8 in the matrix A 
(parameters X ar , with a = b, f and r = 1, 2, and their complex conjugates). 

For the assumed ordering of the eigenvalues e ar , the range of the absolute 
values of the ordinary parts of the commuting polar variables is specified by 
the inequalities 

< ord| Af,i| < ord|Afe 2 | < oo , 
< ord|A/i| < ord|A /2 | < 1 , 

< ord|x pQ ,| < 1 . (49) 

The phase angles of the ordinary parts are allowed to take on all values 
between and 2ir. The integration measure d/i([/A[/ _1 ) in polar coordinates 
is obtained from the Berezinian We get 



dniUAU- 1 ) = d/x(l7)d//(A) . (50) 

Here d/j,(U) denotes the measure for integration over the matrices U, a prod- 
uct of measures for integration over the matrices U x p = U x Up and C/ 7 , 

dfi{U) = dtx{U xf3 )dfi{U y ) . (51) 



14 



We have defined 

muxp) = n ^ pa dx; a • n d(3 pq d(3; q ■ na - ^w^m , (^) 

pa pq p 

dfi(U 1 ) = l[d lpr d 1 ; r . (53) 

pr 

The quantity d/i(A) denotes the measure for integration over the matrices A, 

d/i(A) = ndA ar dA^ r • n(Kii 2 - i^i 2 ) 2 • nd^i 2 + 1 vit 2 • (54) 

ar a rr f 

Eq. ([54]) shows that the measure d//(C/AC/ _1 ) contains a nonintegrable singu- 
larity. It is located at ord( | A^i 1 2 + 1 Ayi | 2 ) = 0, i.e., at the surface of the domain 
of integration. This singularity causes the occurrence of an additional term, 
the Efetov-Wegner term. 



2.2 The Efetov-Wegner Term 

For pedagogical reasons, we first exemplify the treatment of the singularity 
for a simpler case. Indeed, the same problem arises in the integration over the 
four-dimensional matrices T a belonging to the coset space U(l, 1/2) /U(l/1) x 
U(l/1). We return to the full problem below. 

The properties of the four-dimensional matrices T a are the same as those 
of their eight-dimensional counterparts T. The coset space consists of the 
matrices T a satisfying T" 1 = s a T^s a , now with s a = diag(l, 1, —1, 1), modulo 
the subgroup of matrices which commute with L a = diag(l, 1,-1,-1). (We 
label the rows of graded matrices of dimensions 4 and 2 by the indices (p, a) 
and a, respectively.) The matrices have the form 

rp _ | + tfpfl ^12 

a \ t a 2i v^T+^P 
where the two-dimensional blocks an d ^21 are related to each other by 

*2i = U^t (56) 

with k a = diag(l, — 1). The matrix elements of these blocks contain the 8 
Cartesian coordinates of T a . Again, the range of the commuting coordinates 
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is defined by the condition that the ordinary part of 1 + £12^21 be positive 
definite. We consider the integral 

J = Jdfi(T a )F(T a ) (57) 

of a function F(T a ) over the coset space with the measure 

d/,(T a ) = l[d[tl 2 ] aal d[t^ . (58) 



Polar coordinates are introduced |TT| in the same way as described above. 
This yields T a = U a K a U~ l . The anticommuting coordinates reside in the 
eigenvector matrix U a , and the commuting ones, in the "eigenvalue matrix" 
A„. We have 



U a = '^ ° , A a = Vl + AaA a *>-\, 



where u pa and A a denote the two-dimensional matrices 



(60) 

and where A a = fc a A*. Moreover, 

t\ 2 = U la X a U2a, *21 = «2aA a «ra • ( 61 ) 

For the range of the ordinary parts of the commuting polar variables, we 
have < ord|A(, a | < 00, < ord|A/ a | < 1. The integration measure in polar 
coordinates is given by 

dfiiUaAaU- 1 ) = d^(U a )dfl(A a ) , 



p 

d//(A a ) = IJdA aa dA* a • (|A 6a | 2 + |A /a | 2 )- 2 . (62) 

a 

The last line in Eq. (|62]) shows that the measure again contains a nonin- 
tegrable singularity located at ord(| At, a | 2 + |A/ a | 2 ) = 0. Thus, the integral 
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cannot be done by straightforwardly substituting T a = U a k a U~ l in F and 
integrating with the measure (0). This can be seen by the following ex- 
ample. Suppose that F depends only on the eigenvalues Ab a ,A/ a , and does 
not vanish at |Aft a | 2 + |A/ a | 2 =0. Despite the fact that the integral over F 
in the Cartesian coordinates is well defined, the straightforward integration 
in polar coordinates gives an indefinite expression of the type ■ oo, with 
and oo resulting from the integration over the anticommuting and com- 
muting polar variables, respectively. To derive the correct formula for the 
integration in polar coordinates, we proceed in the following way. We start 
with the Cartesian coordinates, exclude from the domain of integration an 
infinitesimal neighbourhood of the singularity of d/i(?7 a A a ?7~ 1 ), and express 
the integral X as the limit 

X = lim fdfi(T a )F(T a )6(v(T a ) - e) . (63) 

Here e is real and positive, and v(T a ) represents the ordinary part of |Ab a | 2 + 
|A/ a | 2 , as given in terms of Cartesian coordinates by 

v(T a ) = [t a 12 ] bb [t a 2iU ~ [tUfMxhf ■ (64) 

Since the integral on the right-hand-side of Eq. ( |63"D is over the region where 
the measure &n{U a k a U~ l ) is regular, this integral can be evaluated by the 
straightforward change to polar coordinates. Thus we have 

J = lim [dniUAaU-^FiUAaU-yiviUAaU- 1 ) - e) . (65) 

We write the function v = v c + v n as the sum of the part v c = \Xb a \ 2 + |A/ a | 2 
which contains only the commuting variables Xb a and A/ a , and the nilpotent 
complement v n , 

Vn = (|A fea | 2 - |A /a | 2 )(7i a 7i a - 72a72a) 

+2*(A ba A} a 7 2a7 * a - A /a A* a7la7 *J 

-2 (|A ba | 2 + |A /a | 2 ) 7 la7l*a72a72a • (66) 

We expand 9(v c + v n — e) in a Taylor series at v c — e. The integral fl65|) 
decomposes into two terms, 

J = l v + l w . (67) 
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The volume term has the form 

ly = lim JdfJiiUaAaU-^FiUaAaU-^eiVc - s) , (68) 

and the boundary term or Efetov-Wegner term is given by 
l w = lim Jdii{U a k a U- x )F{U a k a U- x ){5{v c - e)v n + U\v c - e)v 2 n ) . (69) 

Terms of higher order in v n vanish. To calculate X\y, we introduce new 
commuting coordinates p, (fi, (fib and 0/ by 

A 6a = Vpcos0e^, A /a = y/pwyQj+t. (70) 

Since v n is linear in p, the singularity p~ x in the measure is cancelled by the 
factors p and p 2 arising in the terms due to the Taylor expansion in Eq. (|69|) . 
Because of the terms 5{p — e) and p5'(p — e) in the integrand, the function 
FiUaKJJ' 1 ) can be set equal to its value at T = 1. With 

v 2 n = -2p 2 llal \ al2al * a) Jdp(U a )vi = -^pi, j = 1,2, (71) 

we find that 

l w = F(l) -^— 2 lim / dpdcfidfodfa sin 20 8(p - e) = F(l) . (72) 

In summary, we have 

J = lim Jdp{U a KU~ l )F{U a KU- l )6{\\ ba \ 2 + |A /a | 2 - e) + F(l) . (73) 

This is the desired result. It expresses the integral as a volume term and a 
surface contribution. 

The limit in the volume term can be taken after using the Taylor expan- 
sion of the function FiJJaKaU^ 1 ) in the Grassmann variables 7 pa ,7p a - Except 
for the term of zeroth order, F^UaKaJJ' 1 ), all higher terms in the expansion 
tend to zero as A a tends to the unit matrix. We assume that these terms 
vanish as (|Afe a | 2 + |A/ a | 2 ) a with a > 1. Then, only the expansion term of 
fourth order F^UaAaU' 1 ) contributes. Working out the limit yields 

J = Xy + 2 W = JdpiUaAaU-^F^UaAaU- 1 ) + F (l) . (74) 
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We apply this method to an integral over a function F(T) depending on 
the eight-dimensional matrix T, 

1 = Jdfi(T)F(T). (75) 

From Subsection 2.1 it follows that the matrix T is given by the product 

T = U x(3 dmg(T 1 ,T 2 )U-p 1 , (76) 

where 

yl + ^12^21 ^12 I „ _ •]_ 2 (77) 



\ ^21 VI + ^21^12 / 

are the elements of the coset space U(l, l/2)/U(l/l) x U(l/1) with polar 
coordinates 7 pr ,7pr. and A ar ,A* r (cf. Eqs. ( |j9|) and flBTf )), 



T r = U r A r U r -\ U r = I " lr 1 , A 



\ . / A/i+ArAr A r 



M 2r / ' ' \ A r \/l + A r A r 



»ir = e*P^ j, « ap = «p^. 7 . J, A P =^ Q A/ J -(78) 

In Eq. (|76|) , we now regard the matrices T r as functions of t\ 2 and ^i- This 
defines a new parametrization of T in terms of the polar coordinates of U x p 
and the Cartesian coordinates of Ti and T 2 , with integration measure 

d(U x p diag(T 1; T 2 ) U-£) = d/i(^ /3 )d/i(T 1 )d /U (T 2 )m(A). (79) 

Here &ii{U x p) denotes the measure fl5"2"|) for integration over the matrices C/^, 
d/i(T r ) denotes the Cartesian measures for integration over the matrices T r 
(cf. Eq. ©), 

d/i(T r ) = n d K 2 ]-' d K 2 ]^ , (80) 

and m(A) denotes the factor 

m(A) = n(l A -i| 2 - l^| 2 ) 2 • II (M + |A /r f )- 2 . (81) 

This measure does not contain any nonintegrable singularity in the domain 
of integration. Hence, we can write 

X = Jdfi(U xP diag(7i, T 2 ) U- p l )F{U xP diag(T 1; T 2 ) C^ 1 ) . (82) 
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The source term fl35| ) vanishes at T = 1 where Q = L. Thus in the case 
of Eqs. ( p4|) and (|89D , the function F{U x p diag(Xi,T 2 ) U~p) tends to zero, 
F(l) = 0, when T\ and T 2 tend to the unit matrix. In the integrals over T\ 
and T 2 , we change to polar coordinates, using the rule ([74]). This yields 

X = Xy + X{y , (83) 

X v = Jd^UAU'^F^UAU- 1 ), 

l w = Jdfi(U xP )dfi(U 2 A 2 U^ 1 )F 04 (U x pdia,g(l,U 2 A 2 U^ 1 )U- p 1 ) . (84) 

Here F nin2 (U AU~ l ) denotes the part of F{U AU~ l ) which is of order ri\ in 
7 P i,7pi and of order n 2 in 7 P 2,7p2> an d d/i(?72A2^2~ 1 ) denotes the integration 
measure in polar coordinates for the matrices T 2 , 

d/i(f/ 2 A 2 f/ 2 - 1 ) = dfi(U 2 )dfx(A 2 ), 

MU2) = I[d% 2 d-f* 2 , 

p 

d/i(A 2 ) = \[dK 2 d\* a2 ■ (\\ b2 \ 2 + |A /2 | 2 )- 2 . (85) 

a 

This is the final result. 

We simplify the notation by writing the integrals Xy and Tw as 

X n = / d^i{UAU- l )F{UAU- 1 ) , Q = V, W . (86) 
Jn 

For Q = W, the matrices C/, A and the measure d/i(C/AL r_1 ) are by definition 
given by 

U = U xP diag(l, U 2 ) , A = diag(l, A 2 ) , 

dfiiUAU- 1 ) = dfi(U xP ) d/i(f/ 2 A2^ 1 ) • (87) 

It is always understood that the only non-vanishing contribution arises from 
that part of F which is of the highest possible order in 7 pr ,7p r , i.e., of 8th 
order for Q = V, and of 4th order for Q = W. 
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3 Asymptotic Expansion for Small t 



The behaviour of the correlator C(t, M) = Sg^dg® at small t is described 
by the leading terms of the asymptotic expansion 



C(t, M) = gWgW - gW gW = £ c(n, M) t n - (M/2) 2 , 

n=0 

generated by expanding the exponential in the integrand of Eq. ([34]) in powers 
of t. The expansion coefficients c(n, M) are given by the graded integrals 

c(n,M) =^ n Jdfx(T)D(M,Q)R(Q)S n (Q) , Q = T~ 1 LT, (89) 

where 

AT n = (-) n /(2"n!) , D(M, Q) = detg" M (l + QL) , 

5(g)=trg[(gr 3 ) 2 ], (90) 

and where R(Q) denotes the source term (p5|) . We note that by definition, 
the channel number M assumes even values only. For given M, only the first 
M/2+1 terms of the expansion (^) are meaningful: for higher values of n, the 
factors S n (Q) cause the integrals Q8"5D to diverge (diverging integrals over the 
bosonic eigenvalues). We calculate the first three terms of the expansion. The 
expansion coefficients (|S9| ) can be worked out analytically. The calculation 
relies on symbolic manipulation utilities (we have used those of Mathematica) 



in an essential way and is described in Subsection 3.1. The results are given 



in Subsection 3.2 



3.1 The Coefficients c(n, M) for n < 2 

We change to polar coordinates, use the integration formulae ( |S3D and 
and write 

c(n,M) = Y,cn(n,M) , (91) 

Q 

where 

c n (n,M) =M n f d^UAU- 1 ) D(M,Q) R{Q) S n {Q) , 



n 



Q = UA^LAU- 1 . (92) 
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Both the volume part cy(n, M) and the boundary part cw{n, M) are evalu- 
ated in the same way. We only sketch the main steps. 

Following VWZ, we express the complex variables X ar in terms of absolute 
values and phase angles, 

I 

X ar = i sin -9 ar e^ ar , 9 br = iiV , 9f r = i?/ r , (93) 

with d ar and <fr ar real. Substituting this expression in T = UAU~ l and 
Q = UA^LAU- 1 yields 

Q = UA^LAU^ 1 , A- l LA = ( COs9 lsin ° ) , (95) 

y — z sin # — cos 9 J 

where 

E/ = diag(£ti,-u 2 ) , ui = , u 2 = u 2 . (96) 

The matrices £/, A differ from their counterparts U, A only in their dependence 
on the phases 4> ar . We express the matrices T and Q in terms of the modified 
matrices U = diag(fti, u 2 ) and A, and omit the tilde. It is advantageous to 
parametrize the modified "eigenvalue" matrix A in terms of 

Par = COS 2 -9 ar . (97) 

Explicit expressions for the modified matrices U and A and for the associated 



integration measures are provided in Appendix |TT 
Inverting the matrix 1 + A _1 LAL, we obtain 

( 1 + QL)-^U { 1 + ^LKL)-^ = U\()^ ¥ \^¥y,. (98) 
The term D(M, Q) depends only on the matrix A, 



D(M, Q) = detg~ M (l + A~ X LAL) = D(M, A- 1 LA) , (99) 
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with 



D(M, A _1 LA) 



(/i/i/i /2 ) M (/i 6 i/i 62 )- M « = V , 

for 

/ijfoi" n = w. 



(100) 



Substituting Eq. (0) into the matrices G {rp) {Q) = (1 + QL)- 1 ^ of the 
source term and separating the parts which contain different powers of 

M, we get 

R(Q) = MQ) + MQ) + MQ) , (101) 

where Rf.(Q), k = 2,3,4, denote the source term components 



M 2 

R 2 (Q) = — tig 



r(2). 



r(2). 



r(l), 



M 3 

^ 3 (Q) = -^trg 



x/ (1) («i)x^ (1) («2)x/ (2) K)x^ (2) K) 

+Xl {1) (u 1 )xI {2 \u 2 )xI {2) (u 1 )xI {1 \u 2 y 



M 4 r i r n 

i24(Q) = tt trg X / (1) («i)x/ (1) (« 2 ) trg x/ (2) («i)x/ (2) («2) -(102) 

54 L J L J 



The dependence of Rk(Q) on U is contained in the graded matrices I^ r '\ 



u. 



Vli 



I {r \u p )=u p 



^I^Up , = diag(l, 0, -1,0), J (2) = diag(0, 1, 0, -1) 
their dependence on A in the matrix 

X = tan -0 . 

Substituting Eq. ( j95|) in the coupling term S'(Q) yields 

S(Q)=S n (Q)+S 2i (Q) + 2S i3 (Q), 
where S pp /(Q) are the graded traces 



(103) 
(104) 

(105) 
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Sn(Q) = trg cos9t 3 (ui) cos0t 3 (ui) 
S22(Q) = trg cos 6 t 3 (u 2 ) cos6t 3 (u 2 ) 
SuiQ) =trg sm6»r 3 (iii) sin 6>r 3 (w 2 ) 

which depend on U via the graded matrices 

r 3 (u p ) = u~ l TzU p , r 3 = diag(l, -1, 1, -1) 



(106) 



(107) 



Writing the measure d//(C/AC/ _1 ) as the product of measures for integra- 
tion over U and over A, using Eq. (|101|) and performing first the integration 
over U gives the coefficients Cn(n, M) as the sum of contributions of the three 
source components, 

c n (n,M) =$>[?(n,M) , 



:io8) 



where 



with 



(n,M) = I dfi(A)D(M, A^LA)^ (n, M, A) , (109) 
Jn 

F$\n, M, A) = Af n ( d l i{U)R k {Q)S n {Q) . (110) 
Jn 

The measure d/j,(U) is the product of measures for integration over the matri- 
ces u\ and it2, the source components Rk{Q) and the traces S PP >(Q) depend 
on Up only in terms of the matrices I^ r '(u p ) and ^(tip), and this dependence 



has the simple form shown of Eqs. ( |102j ) and ( |106| ). For all these reasons, 
the evaluation of Fq(ji, M, A) reduces to the calculation of the integrals 

Kq (h,h,i2,h\h,h, ■ ■ -ikmJm) 

= ^ d ^K)[ /(1) K)]iiii[ /(2) K)]i 2 j 2 [ r 3K)]fc 1 h • • • [r 3 (u p )} kmlrn , (in) 

where p = 1,2. From the mathematical point of view, these integrals are 
linear forms of integrals of products of matrix elements of the coset matrices 
Up, u^ 1 over the coset Z7(2/2)/Z7(l) x (7(1) x U{1) x U(l), and could therefore 
be evaluated by a generating function approach similar to that worked out 



for the same integrals over the graded unitary group by Guhr [18, [19 



The integrals Kq and Kq referring to the same set of indices ii,ji,i2, J2, 
kx, lx, . . . , k m , l m are closely related to each other. Starting from the integral 
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Ky\i\, ji,i 2 , h, ■ ■ ■ , k m , lm), making the change to the primed anticom- 
muting variables introduced by f3 lq = i(3' lq , (3{ q = i/3' lq , j lr = ij' lr , 7^ = i^f[ r , 

and comparing the result with the integral Ky\ii,j 1 , «2, J2 1 ^i , h, ■ ■ ■ , k m , l m ), 
we find that the two integrals in fact differ only by a factor, the integral over 
the phases <p ar which is present in Ky^ and absent in Ky\ Apart from an 
additional difference in sign due to the Berezinian, the same is true for the 
integrals K$ and Kyy ■ This simplifies the calculation in an essential way. 
The integration over the anticommuting variables yields the part of highest 
possible order of 

(1 - 2/3 pl ^ 1 /5 p2 /3; 2 )[/( 1 )( Mp )] ilJ J/( 2 )( Mp )] i2j2 [r 3 K)] fclil . . . [rs(u p )] kmlm (112) 

(of 16th order for = V, of 12th order for il = W). After setting x pa = 
sin Cpa eiripa ) the integration over the commuting variables appearing in u p 
simplifies to integrals over phase angles, or over products of powers of basic 
trigonometric functions. Multiplying the products of complementary pairs of 
integrals K^\i\,j\, «2, J2 1 ^i , h, ■ ■ ■ , k m , l m ) by the A-dependent factors stem- 
ming from x an d cos6> and collecting the contributions yields the desired 
F^ h \n, M, A). The result has the form 

4 fc) (n, M, A) = N nk M k ]7>6r - H fr >) ■ U^r ■ Py \n, A) , (113) 

rr ' ar 

F$\n, M, A) = N nk M k {^ b2 - /x /2 ) • Pw(n, A) , (114) 

a 

where P^\n,A) are polynomials in fi ar , and where 

Afnk = (-T(l + 5 k2 )/(2 n+2 n\). (115) 

The polynomials Py^ (n, A) are symmetrical with respect to the interchange 
of /Xfti and /z^, and of /ifi and /i/ 2 . Taking into account this property, we can 
write the final integrals over A as 

c { y\n, M) 

= M nk M k f d/x(A)n^-X", M " 2 - 11(^-^0- Pp\n,A) 

JV r rr' 

]_ roo roo rO rO 

= Af nk M k -J d/i 61 J dfx b2 J d/x/i J d/i /2 x 
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c^(n, M) = A/; fc M fc / d/i(A)/ig- 2 /Ub - 2 M " 2 (^2 - Ji/a) P«(n, A) 

= M nk M k d/W d/i /2 4^ P$(n,A). (116) 

The integrals can be done analytically using the approach described in Ap- 
pendix |6.2| . Collecting the contributions of all three source components finally 
yields the coefficients Ca(n,M). 

The calculation of the volume integrals Cy\n,M) can be simplified by 
using modified polar coordinates where 

T = UAU- 1 , U = diag(«i, u 2 ) , 

Ui = U^UxpUxxUtf, , u 2 = u 2l u w u 2x , (117) 

with u pl ,u p p,u px and having the same form as in the old parametriza- 
tion. The integration measures for both parametrizations coincide. Since 
the matrices u P1 commute with the matrix T3, the matrices t^{u p ) are now 
independent of 7 pr ,7p r - m the integrals Ky , the only part of the integrand 
which then depends on these anticommuting variables is the product of ma- 
trix elements of the matrices I^ r '(u p ), and only the part of highest order 
in 7 P r,7p r contributes. Straightforward calculation shows that this part is 
obtained by replacing the matrices u^I^'u^ in I^(u p ) by the matrices 
(— ) p 27 pr 7* r J (r ), where the J^'s denote the projectors = diag(l, 0, 1, 0) 
and J( 2 ' = diag(0, 1, 0, 1). This substitution simplifies the calculation very 
much. Unfortunately, no similar coordinate transformation has been found 
to simplify the boundary integrals. 



3.2 Leading Terms 

We begin with the expansion term of zeroth order c(0, M). The integrals 
Kq which enter the calculation do not contain any of the matrix elements 
of t 3 (u p ). Integrating over u v we find that Ky Ji, Ja|) = for all 
*i)Ji)*2)j2- Thus Cy (0, M) = for any k, and the coefficient c(0, M) is 
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entirely determined by the boundary terms c^)(0,M). We obtain, after 
combining the contributions from the integrals K^(ii, ji, 12, J2I), that 

P^(0,A) = -fi b2 + H f 2 , 

P^ 4) (0,A) = P^(0,A). (118) 
The integration over A then yields 



f (0,M) = 2(il ^ 2 _ i) =-c^(0,M), 4 ) (0,M) = 4(]i ^ 4 - . (119) 



(2) (3) 

The contributions c^/(0, M) and c^/ (0, M) cancel each other, and the coeffi- 
cient c(0, M) is given by the boundary contribution of the source component 
R 4 (Q), 

c(0, M) = c$(0, M) - \M* = 4(i ^ 2 _ 1} • (120) 

By definition, c(0, M) = C(0, M) gives the variance Var[g(P)] of the dimen- 
sionless conductance g{B). The value (|120|) of Var[g(B)] agrees with the one 
derived by Baranger and Mello using a random S-matrix approach [fZU . 



We turn to the coefficient of the term linear in t. Here, we meet a similar 
situation: All integrals Ky\ii,ji, i 2 , J2I&1, ^1) vanish. With 

P$(l,A) = (8/3) x 

{fi b2 {2 - fi b2 - 3/4) - /4( 2 M&2 - 3)/i /2 - /4(3 - Afi b2 )fi 2 f2 } 

+ (b2 <- /2) , 
P^(l,A) = -(16/3) x (/i fe2 -/i /2 ) 

{1 - /i 62 (l + 3/i 62 ) + 2/ife 2 /i /2 + ^2^/2 + ( b2 ^ f 2 )} > 
P^ 4) (l,A) = (8/3) x 

{/i fe2 (2 - fi b2 - 3/4?) - fi b2 (6 - 3/^2 - 10/i^ 2 )/i /2 - /4( 8 /^2 - 3)//j 2 } 

+ (62 <- /2) , (121) 

the boundary contributions of the first two source components compensate 
each other, 

clJ } (l,M) = - (M 8 2 ^ 1)2 = -c$(l,M), (122) 
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and the coefficient has the value 

c(l, M) = cg?(l, M) = - (M 4 2 ^ 1)2 • (123) 

We finally turn to the term of second order in t. Here, both cy(2, M) and 
cvi/(2, M) contribute because Ify^^i, ji, i 2 , J2I&1, h, k 2 , h) does not vanish. 
The polynomials P#°(2,A) and P$(2,A) are given by 

(2, A) = (128/9) x 

{ [ 8j2 b ijj, b2 (1 - 2(// w + + 4/i W /ife2 + A*ii + a4 ) 

+8(a4 + a4) - 16 (/ i w + a4) + 12 (/4 + /4) ] A*/iA*/2 

-[ftbl^i 7(/ifel + /ife2) - 4(/i 6 i/i 62 + /4 + /4) ) 
-4(^61 + a4) ] /W/2(^/l + 

-4/i W /i b2 [2(/i W + /i 62 ) - 3/i W /i b2 ] (a*/i + M/ 2 ) + 36 a4a4 A*/iA*/2 
+(M<-/1, 62<-/2)}, 
P^ } (2, A) = -(64/45)(^2 - /i/ 2 ) x 

{ 7 - 52^2 + 70/4 - 100/4 + 90/4 

-2/i b2 ( 15 - 70fi b2 + 70/4 - 42/i^ ) /"/a 
-/4( 225 - 404/i b2 + 162/4 ) ^ - 86/4 /ij 2 + (62 <- /2) } , 

Pf (2, A) = (512/9) x 

{ [ 2fl bl fl b2 (fl 2 bl + + yUfei/ifes) + 2(/4 + /4) 

-4(a*M + a4) + 3 (/4l + A4 )] 
+ [A i f>l/ i f>2( /4 + /4 - 8/if,i/if,2 ) + Atfoi + a4] /i/i/// 2 (/i/l + /i/2) 
+/ibl/i& 2 [ 2 (/!{>! + /i 62 ) - 3/i W /i 62 ] (/i/i + A*/ 2 ) 

-(&1<-/1, 62^/2)}, 
P l g ) (2,A) = (256/45) x 

{ /i& 2 (7 - 26/i b2 + 35/4 - 50/4 + 45/4) 

2 



+/i b2 (5 - 35/i 62 + 10// 62 + 25^2 - 42^)^/2 
+A*6 2 (90 - 115//62 + 20/4 - 9/4) A^?2 
+5/4(5 + 3/4,2) At/2 + (62 <-> /2) } , 
p( 4) (2,A) = pf(2,A), 
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P^(2,A) = -(2/45)(// 62 - M/2 )x 

{ 7 - 52fi b2 + 70/4 - 100/4 + 90/4 
+2/x fe2 (75 - 230/x b2 + 230/4 - 138/4) /i /2 
+/4(435 - 676/i b2 + 198/4) /^ 2 

+154/4/4 2 + (62^/2)}. (124) 
The contributions of the three source components to cn(2,M) are 

2Af(M 6 - 9M 4 + 18M 2 + 18) , 
4AT(4M 6 - 12M 4 + 45M 2 - 9) , 

27VM 2 (11M 2 - 39) , 
-27VM 2 (9M 4 - 22M 2 + 69) , 



c^(2,M) = 
c$(2,M) = 

4 3) (2,M) = 
c$(2,M) = 



c^(2,M) = J\fM 2 (M 6 - 9M 4 + 18M 2 + 18) , 

c$(2, M) = -7VM 2 (M 6 - 18M 4 + 51M 2 - 90) , (125) 

where 

M = 16/(3(Af 2 - 1) 2 (M 2 - 4)(M 2 - 9)) . (126) 

Similarly as in the terms of zeroth and first order, the total contribution of the 
source components -R 2 (Q) and R?,(Q) is equal to zero, c^(2, M)+c^(2, M) = 
0, and we find that 

u), N us, . 16M 2 (3M 4 - 11M 2 + 36) , x 
c(2, M) = c<?(2, M) + eg>(2, M) = ^ _ L_ _ _ _ L . ( 12 7) 

Collecting the terms of all three orders, we get the leading part CAE(t,M) 
of the asymptotic expansion (|8~8"D in the form 

M 2 4M 3 
C ^ (t ' M) " 4(M 2 - 1) " (M 2 - l)^ * 



16M 2 (3M 4 - 11M 2 + 36) 2 



(M 2 - 1) 2 (M 2 -4)(M 2 -9) 



t 2 . (128) 



(For M — 2, only the first two terms are meaningful.) This is our final 
analytic result. 
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3.3 Comparison with the Lorentzian form 



Efetov || has shown that in the limit M> 1, the autocorrelation function 
has the form of the square of a Lorentzian, 



C SL (t,M) = C(0,M) 



1 + (t/t d ) 



(129) 



where td denotes a decay (shape) parameter. Thus, it is of interest to com- 
pare our analytical result (|128|) with the Taylor expansion of the squared 
Lorentzian at t — 0. It is given by 



C SL (t,M) 



M 2 



4(M 2 



4» 



+ 3(f) 2 + 



(130) 



For M 3> 1, the two linear terms become identical if we put td = M/8, and 
the quadratic terms then coincide in this limit. This confirms the expected 
agreement for M 1. However, for the small channel numbers M typical 
for the experiments mentioned in the Introduction, no choice of td exists 
for wich the linear and quadratic terms in the expansion of the squared 
Lorentzian would agree with the form (|128|) . Put differently, with td = M/8, 
the decrease of the squared Lorentzian near t = is much smaller than that 
of the autocorrelation function ([128]) . A numerical test for the asymptotic 
expansion (|128|) is presented in the next Section. 



4 Numerical Results 

To allow for a comparison between the squared Lorentzian and the auto- 
correlation function of our random matrix model also for larger values of t 
but small M, we have calculated the latter function numerically. Here we 
present only the main results of this study. A more detailed presentation of 
the numerical method and the results will be published elsewhere. 

After rescaling and the introduction of the now more convenient param- 
eter b = k(B^ — 5 (2) )A/0 O , b 2 = At , the quantity of interest takes the 
form (cf. Subsection |L~ll ) 

2 (131) 
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The conductances g^ r ' are given in terms of the S matrices S^ r \ The latter 
have the form 



S (r) = 1 - 2iitW[D {r) }~ l W T , D {r) = -H {r) + iW T W , 



W Cil = 5 C „ , H$H$, = , C = I,II. (132) 

The ensemble average was performed by repeated drawings of the random 
matrices H^ 11 ^ from a random number generator. The results for M = 
2,4 and 10, calculated with matrices of dimension N = 50, are shown in 
Fig. 1. The solid line shows the best fit by a squared Lorentzian written in 
the form 

r i i 2 

C SL (b,M) = C(0,M) - - -j— , (133) 

with (7(0, M) = M 2 /[4(M 2 - 1)]. The values of the shape (decay) parameter 
hd are presented in Table 1. For comparison, the table gives also the values 



of b d suggested by the results of Efetov || and Frahm jl(J 



The table shows that the best values for bd are very close to Efetov's value 



M/2. For all values of M considered, the autocorrelation function (7(6, M) 
agrees rather well with the squared Lorentzian. As expected, the agreement 
improves with increasing M. For M — 2, the deviations of (7(6, M) from 
the best squared-Lorentzian fit Csh{b, M) (with bd = 1.1) are shown in more 
detail in Figs. 2 and 3. Even in this case, the deviations do not exceed 5% in 
magnitude. For M > 2, the best fit curve for CsL{b, M) lies within the error 
bars of the numerical calculation. 

As a test of the asymptotic expansion ( |128| ), we have calculated numer- 
ically also the second derivative d 2 C(b, M)/db 2 at b = for the two lowest 
M values, M = 2 and M — 4. We assumed that taking the derivatives can 
be interchanged with taking the average. The derivatives of the integrand 
in Eq. ( |131| ) were analytically done by applying to the resolvents the 
formula 

d n 1 1 / 1 \ n 

dx n A + xB = n!(_)n A + xB \ B A + xB ) ' (134) 
with A, B denoting matrices independent of x. Table 2 compares the second 
derivatives with the values derived from the asymptotic expansion, and with 
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the corresponding derivatives of CsL(b, M). Expressed in terms of b, the 
expansion takes the form 



^ b,M ^ 4(M 2 -1) (M 2 -l) 




v 1 b ± 

(M 2 - 1) 2 (M 2 -4)(M 2 -9) 



M 2 (3M 4 - HM 2 + 36) 



(135) 



The values derived from the asymptotic expansion agree very well with the 
numerical results for both M = 2 and M — 4. This is not true for the 
derivative of the squared Lorentzian which for M = 2 differs considerably 
from the numerical result. 

The same calculation of the fourth derivative <9 4 C(0, 2)/db 4 yielded a di- 
vergent result. This fact caused us to analyse the second derivative d 2 C(b, 2) / 
db 2 at small |6| in greater detail. The result is shown in Fig. 4. In contrast to 
the second derivative of the squared Lorentzian, which at the origin rises with 
the second power of b, the second derivative of C(b, 2) seems to rise with the 
absolute value of b. The derivatives of this term linear in |6| then generates 
a 5-function singularity in the fourth derivative. In the power spectrum of 
C(b, 2), i.e. in the Fourier transform, such a singularity manifests itself in an 
algebraic decay for large k. 

In summary we have shown that at the 5% level of accuracy, there is no 
difference between our results and a squared Lorentzian even for the smallest 
M value, M — 2. Closer inspection of the autocorrelation function for M = 
2 and at the point b = shows, however, that there is strong evidence 
for a singularity of the fourth derivative, caused by non-analytic behavior. 
This is an interesting and unexpected result for which we have no physical 
explanation at present. 



We have investigated the magnetoconductance autocorrelation function for 
ballistic electron transport through microstructures having the form of a 
classically chaotic billiard. The structures were assumed to be connected to 
ideal leads carrying few channels. Assuming ideal coupling between leads and 
billiard, we have described this system in terms of a random matrix model. 



5 Summary and Conclusions 
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The autocorrelation function depends only on the field parameter t, spec- 
ified by the field (flux) difference, and on the channel number M. Using 
the supersymmetry technique, we have analytically calculated the leading 
terms of the asymptotic expansion of the correlation function at small t. 
To this end, we have used integral theorems obtained by applying Berezin's 
method of boundary functions. Using a generalization of the standard po- 
lar coordinates to parametrize the coset space, we succeeded in identifying 
and evaluating both volume and boundary (or Efetov-Wegner) terms. We 
believe that the method developed in this paper is of general interest for the 
supersymmetry technique, and we hope that it will be helpful in other cases. 

We have shown that the first two terms of the asymptotic expansion of 
the autocorrelation function are entirely given by the Efetov-Wegner terms. 
This result is likely to be of general importance: Given some seemingly very 
natural choice of integration variables, the boundary terms may easily yield 
a or the major contribution. 

For large M, the asymptotic expansion agrees with the corresponding 
small t expansion of the squared Lorentzian suggested by semiclassical theory. 
For small M, differences exist. These were studied further by combining our 
analytical work with numerical simulations. For the smallest value of M, 
M = 2, the difference between the autocorrelation function and the best 
squared-Lorentzian fit exists but does not exceed 5% in magnitude. This may 
seem irrelevant. However, a study of the derivatives of the autocorrelation 
function yielded further evidence for a statement suggested by the analytical 
form of the asymptotic expansion: The autocorrelation function seems to be 
non-analytic in t at t — 0. We find this result surprising. It suggests the 
occurrence of non-analytic behavior also in other correlation functions where 
the consequences may even be observable. 
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6 Appendices 



6.1 Matrices and Measures 

In the volume integrals Ty, the matrices U nad A have the form 



U 



ui 
u 2 



\ — ( cos \® * s * n f ^ 
I % sin |0 cos \d 



(136) 



where # denotes the diagonal matrix 6 = diag(ii?{,i, ^62, 1^/2), an d where 
■Ui and w 2 denote the matrices 



Ml = U^UipU^Ufr U 2 = U 2x U 2 pU 2l , 

obtained by multiplying the factors (we set x pa = sm( pa e" lpa ) 



(137) 



u 



px " 



u 2 p-- 



u 2l - 



( cos Cpb sin C P bZ ir,p " 

- sin Cpbe" 1 ^ 6 cos ( pb 



\0 

fl + lPnPt, 






/ 1 — \P 2 iP 2 \ 



V^2*l 

/l + i 7 n7n 



1 + \M 






1 

t(3. 




22 



22 



\ 



cos C P f sin C P /e l,?p/ 

- sin Cp/e - **' cosCp/ 

Ai \ 

A 2 

l-|/3i2^ 2 
1 - i/Jn^ ; 

i/3 21 \ 
i/3 22 

1 + 1/322/9^2 











7fi 
VO 



*721 

VO 



1 + |7l27?2 



7*2 



7n 




1 - 0-722722 



^722 



712 

1-I7n7n 
1 - 1-71271*2 J 

il2i 







v 

1 + |72i72i 
1 



^722 



(138) 



1 + i722722 J 
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and Ufj, = e^, with = diag(0&i, 062, <f>fi, 0/2)- The corresponding integration 
measure is 

d^UAU- 1 ) = d/i(C/)d//(A) , (139) 

where d[i(U) denotes the measure for integration over the matrices U given 
by 



a 



p) ' 



d/j,(ui) = d / u(Mi x )d / u(Mi /3 )d/i(Mi 7 )d / u(0), 
d/i(u 2 ) = d^{u 2x )d^{u w )d^{u 2l ) , 
djJi(up X ) = Y[ d( pa dr] pa sin 2( p 



ipa i 



dfj,(u^) = J[ d0 ar , 

ar 

q p 

d/i(w P7 ) = PJd7 pr d7p r , (140) 

r 

and where d/j,(A) denotes the measure for integration over the matrices A, 

d/i(A) = Yl <W ■ YliPod - A^) 2 ■ YiiVbr - ^fr')' 2 , (141) 
ar a rr' 

with /j br = cosh 2 ^$br and \if r — cos 2 |$/r- The domains of integration over 
(the ordinary parts of) ( pa extend from to tt/2, those of integration over 
Vpa, (par from to 2n , and those of integration over ^ and fif r from 1 to oo 
and 1 to 0, respectively, with \i al < ji a2 . However, as discussed in Subsection 



34], calculating the volume integrals can be simplified drastically by using 
a modified parametrization where the matrices u\ and u 2 are given by the 
products 

ui = u^UxpUixUj,, u 2 = u 2 ^u 2/3 u 2x . (142) 

The integration measure in the modified parametrization has the same form 
as in the old parametrization. 

In the boundary integrals Xw, all variables labelled by r = 1 are to be set 
equal to zero: the matrices u px and u p /3 appearing in U have again the form 
shown in Eqs. (|138|) , whereas the matrices u\ 7 , u 2l , and 9 simplify to 
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Win 



«2n 



/I 











\ 





1 + hl27l*2 





712 










1 







\o 







1 - 


|7l27l*2 / 


/i 











\ 





1 - ^722722 





^722 










1 







Vo 


*722 





1 + 


|722722 / 



= diag(0, fe2 , 0, /2 ) , 
= diag(O,n? b2 ,O,tf /2 ) . 



(143) 



In the integration measure dfi(UAU x ), the measures for integration over the 
matrices u pl , and A simplify to 



d/i(w P7 ) = d7 p2 d7p 2 , d^(u^) = JJ d0 o2 , 

a 

d//(A) = JJd^ • (/i62 - /U /2 )" 2 . 



(144) 



6.2 Integration over Eigenvalues 

The integration over A of F$\n, A, M) leads to the two-dimensional integrals 



I(p2,q2) = / d/i 62 / d/j, f 

J 1 «/ 



M-2+g2 1 

1 

2 M+2-p2 ,, 



(145) 



with p2, q2 nonnegative integers. The integrals converge if M + 2 > p2 and 
M + q2 > 1; since M > 2 for physical reasons, the second condition is always 
satisfied. With x = /if 2/ 1 [Xm and y = /i/ 2 , we can write 



l( P 2,q2)= / dx- / dyy 

JO 1 — X Jx 



.p2+<j2-4 



(146) 



For p2 + q2 7^ 3, the integration over y gives 

1 n 



I(p2,q2) 



dx 



x M-2+q2 _ x M+l-p2 



X 



3-p2-q2Jo 
3 _ p 2-g2 (4>(M + 2 — p2) — tp(M — 1 + g2)) , (147) 



36 



where ip(z) denotes Euler's if) function pi| . Making use of 

1 

k 



^(p + i) = Er + #) 

fc=i 



(148) 



yields 



I(p2,q2) 



I(p2,q2) 



M+l-p2 1 

y - 

k=M-l+q2 ^ 
M-2+q2 * 

p2 + g2-3 k= ]t2- P 2 * 



3 - p2 - q2 
1 



for p2 + q2 < 3 
for p2 + q2 > 3 



For p2 + g2 = 3, the integration over y leads to 



J(p2,3-p2) = - / 



1 x Af+l- P 2 lna , 



dx 



1 — a; 



C(2,M + 2-p2) 



(149) 



(150) 



where C( z j (?) denotes Riemann's ( function f22|. Substituting the explicit 
expression, we get 



I(p2,3-p2) 



^2 M+l-p2 1 



6 



E 

k=l 



k 2 



(151) 



The integration over A of Fy(n,A,M) leads to the four-dimensional inte- 
grals 



I(pl,p2,ql,q2) 

■, 1 M-2+ql M-2+q2 „ / N9 

[°°A f°°A f A I A ^/l /jfg Ha (^l-/i Q2 ) 2 

- / d/i W / d/i b2 / d/i/i / d/x/2 M+2 _pi M+2- P 2 n r„ n — T" ' 

Jl Jl JO JO /i w T fJ- b2 llrr'U^&r — A*/r'J 

(152) 

with pl,p2, gl, g2 nonnegative integers. Using the identity 



Ha (fJ'al — ^a2) 2 _ Pbl ~ Pb2 Pbl ~ Hb2 l^bl ~ ^b2 Hbl ~ Hb2 
Hrr'i^br ~ fl>fr') Hbl ~ M/l l^bl ~ A*/2 Hb2 ~ /ffl ^62 — A*/2 

(jl b l — p b2 ) 2 (Hbl — Hb2) 2 



(A*61 - A*/l)(A*62 - (A*61 - M/2)(A«62 ~ A*/l) ' 

(153) 
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we can express this integral in terms of the two-dimensional integrals ( |145|) 
just considered, as seen from 



and 



M-2+ql M-2+q2 



f°° f°° f f A*/i A*/2 

/ dfl bl / d^ 62 / d/i/i / d/i/2 M+2 -pi M+2-p2 



/^/l 0/2 f'bl ~ /^62 

M+2-pl M+2-p2~ " 
061 062 0" ~ 0/1 

J(pl + l,gl) 



(M+ 1 - P 2)(M- 1 + g2)" 

1 r^(pl,?l), (154) 



(M-p2)(M- l + g2)" 



/•oo /-oo /-l /-l ^ 

/ d/i W / d/i b2 / d/i/i / d/i/2 - j m+2 
Ji Ji Jo Jo « k1 



Af-2+gl M-2+q2 , 

A*/2 (061 — 0^2 J 



061 ' ~ P Vb2 +2 ^ (0 fcl ~ 0/0 (062 - 0/2) 



= J(pl + 2, gl) J(p2, g2) - 2Z(pl + 1, ql) J(p2 + 1, g2) 
+X(pl,gl)X(p2 + 2,g2) . (155) 
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Table 1: The value of the shape parameter bd for the best squared-Lorentzian 
fit CsL{b, M) for the channel numbers M = 2, 4 and 10. The last two columns 
present the bd values \J M/2 and J (M + l)/2 predicted by Efetov |§ and 
Frahm JTU] respectively. 



M 



M/2 



V / ( M + 1 )/ 2 



2 
4 
10 



1.10 
1.44 
2.20 



1.00 
1.41 
2.24 



1.22 
1.58 
2.35 



Table 2: The value of the second derivative of the correlation function 
C(b,M) at b = for M = 2 and M = 4. The values obtained from the 
asymptotic expansion Cab(&, M) and from the best squared-Lorentzian fit 
Csi{b,M) are shown for comparison. 



M d 2 C(0,M)/db 2 d 2 C AE (0,M)/db 2 d 2 C SL (0,M)/db 

2 -1.76 -1.78 -1.08 

4 -0.52 -0.57 -0.50 
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C(6,M) 



0.4r 

0.35- 




Figure 1: The correlation function C(b, M) as a function of the magnetic-field 
parameter b for the channel numbers M = 2,4 and 10. The solid lines show 
the corresponding best squared-Lorentzian fit Csi{b, M). The calculated 
values of the correlation function are lying in the middle of the error bar 
intervals. 
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6C(b,2) 

O.Olr 



.0075 
0.005 
0.0025 



■0 . 0025 
-0 .005 
■0 . 0075 
-0 . 01 L 





Figure 2: The deviation 5C(b,2) = C(b, 2) - C SL (b,2) of the correlation 
function C(b,2) from its best squared-Lorentzian fit CsL(b,2) as a function 
of b. The gray band connecting the error bars is intended to help the eye. 
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Figure 3: Plot of log 10 C(b, 2) as a function of log 10 b. The solid line refers to 
the best squared-Lorentzian fit, the dashed line to the b~ 2 law, for compari- 
son. 
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Figure 4: Graph of the first two derivatives of C(b, 2) as a function of b. 
The upper curve shows the derivative C'(b, 2) = dC(b, 2)/db, the lower curve 
the function C"(6,2) = <9 2 C(6, 2)/<% 2 scaled by the factor 1/2. The dashed 
and the solid line show their respective counterparts for the best squared- 
Lorentzian fit CsL^b, 2). 
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